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By using tlie Feynman-Hibbs prescription for the evolution amplitude, we quantize the system of a 
damped harmonic oscillator coupled to its time-reversed image, known as Bateman's dual system. 
The time-dependent quantum states of such a system are constructed and discussed entirely in the 
framework of the classical theory. The corresponding geometric (Pancharatnam) phase is calculated 
and found to be directly related to the ground-state energy of the ID linear harmonic oscillator to 
which the 2D system reduces under appropriate constraint. 
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I. INTRODUCTION 

From the very outset of quantum theory a tremendous 
effort has been devoted to answering the question where 
is the boundary between the classical (macroscopic) and 
the quantum (microscopic) world, or in other words, 
when a quantum system starts to behave classically. The 
correspondence principle introduced by Bohr in the end 
of 20 's served as a heuristic prescription to construct 
quantum mechanics. Roughly speaking, the quantum 
theory should approach the classical theory in the limit of 
large quantum numbers. However, in general, this limit 
is quite subtle On the other hand, a typical state- 

ment in majority of standard textbooks is that classical 
mechanics applies in the limit ft — )■ 0, and this paradigm 
has become over the years a starting point for a host 
of semi-classical asymptotic treatments Yet, still 

many precautions should be taken when it comes to non- 
physical infinities at caustics, boundary layers analysis 
or to connection rules . In terms of canonically con- 
jugated variables, say p and q, the connection between 
the classical and quantum descriptions has been estab- 
lished by Ehrenfest's theorems p]J3,p| which give the law 



of motion for the mean values of the operators p and q 
in a form emulating the classical equations of motion. 
However, only for systems with quadratic Hamiltonians 
is this correspondence exact, i.e. the mean values (p) 
and (g) follow classical phase-space trajectories [Q. In 
particular, the system of quantum harmonic oscillators 
with bilinear coupling can be completely characterized 
using classical trajectories because the quantum Wigner 
function satisfies the classical Liouville equation 

In the present article, the first of the series, we attempt 
to shed some further light on the quantum-classical rela- 
tion by studying the quantization of the 2D system of a 
damped harmonic oscillator coupled to its time-reversed 
image, first introduced by Bateman |p^ . Bateman's dual 
system received considerable attention in the past as it 
represents a simple explicit example of a dissipative sys- 
tem which could be tackled by means of canonical quanti- 
zation jri|-p^. However, the Quantum Mechanics (QM) 
of this system is plagued with many conceptual problems 
(e.g., the wave functions cannot be normalized in 
the usual manner, the Hamiltonian is not self-adjoint 
and represents the energy only for a restricted set of dy- 
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namical solutions) and it was shown |Q that a consistent 
quantization can only be achieved in Quantum Field The- 
ory (QFT). More recently, Bateman's system has been 
studied in connection with (Chern-Simons) gauge theo- 
ries , as an example of an exactly decoherent system 
Ip^ and as a toy model for the recent proposal by G.'t 
Hooft about deterministic QM p6|jl7| ]. These aspects as 
well as the QFT of the Bateman system will be the ob- 
ject of future papers. The aim of this paper is to present 
a thorough analysis of the QM of this system which will 
be the basis for the next papers: our results include the 
time-dependent wave functions and the geometric phases 
associated to them. We also study the reduction of the 
2D Bateman system to the ID linear harmonic oscillator 
(l.h.o.). 

An outline of the paper is as follows: 

In Section II we quantize Bateman's system by using the 
Feynman-Hibbs prescription for the time-evolution am- 
plitude (kernel) [lj-|2^: this allows us to avoid the pit- 
falls of canonical quantization. We show that the kernel 
is fully expressible in terms of solutions of the classical 
equations of motion and that it is invariant with respect 
to the choice of the fundamental system of those solu- 
tions. This might be viewed as a two-dimensional exten- 
sion of an existing one-dimensional result An im- 
portant ingredient of the kernel calculation is the fluctua- 
tion factor . This is calculated by employing the Van 
Vleck-Pauli-Morette determinant technique which al- 
lows us to avoid a direct manipulation of the Schrodinger 
equation. Instead, only the phase-space structure of the 
classical solutions (e.g., Lagrangian manifold) is used. 

In Section III we use Mehler's formula [ pT|j22| for a spec- 
tral decomposition of the kernel in order to obtain the 
time-dependent wave functions in hyperbolic radial co- 
ordinates (r, u) . They are expressed in terms of general- 
ized Laguerre polynomials and are shown to satisfy the 
correct time-dependent Schrodinger equation. The ex- 
plicit form of the wave functions uncovers the root of the 
difhculties connected with the canonical quantization - 
the unboundedness of Bateman's system in the hyper- 
bolic angle u. Use of the radial kernel [|l9| allows us 
to factorize away the explicit u dependence. The "ra- 
dial" wave functions then correctly fulfil both orthonor- 
malization and completeness relations. In addition, the 
radial wave functions satisfy the radial, time-dependent 
Schrodinger equation with the Hamiltonian Hi. For the 
"azimuthal" quantum number I = ±i, the latter turns 
out to be formally identical with the Hamiltonian of the 
ID l.h.o.. 

To better understand the structure of the wave functions, 
we focus our attention on the algebraic setting of Bate- 
man's Hamiltonian. Identifying the dynamical group as 
SU (1,1), we are able to pinpoint the structure of the 
ground state, which turns out to be a (squeezed) coherent 
state. The aforementioned peculiar behavior of the wave 



functions is then attributed to the remarkable properties 
of the SU {1,1) group representations: the unbounded- 
ness of Bateman's system in the variable u could be seen 
as a consequence of the non-existence of a unitary irre- 
ducible representation of SU{1, 1) in which the genera- 
tor J2 would have at the same time a real and discrete 
spectrum [ p3[ . The requirement of discreteness of the 
J2 spectrum then leads to an effective non-hermiticity 
and oddness of J2 under time reversal. To accommodate 
this point in the Feynman-Hibbs kernel prescription, a 
new inner product has to be defined. In fact, one of 
the merits of the presented method is that it naturally 
provides the consistent inner product for wave functions 
(the rather artificial and involved method of Racah is 
not needed). In the case when we restrict our attention 
to the stationary quantum states, the connection with 
existing canonical quantization results [pl]-|l3| is readily 
established. 

With the (full) time-dependent wave functions at hand, 
we are able to calculate in Section IV the exact geomet- 
ric (Pancharatnam) phase for Bateman's dual system. 
We find that Pancharatnam's phase is explicitly h inde- 
pendent and consists of three autonomous contributions: 
overall ground-state fluctuations of p and x gathered dur- 
ing the period of evolution and the Morse index. 

We then show that the (full) wave functions become pe- 
riodic in configuration space when the hyperbolic angle u 
solves the classical equations of motion. In this case the 
period of the wave functions matches the inverse of the re- 
duced frequency of the original Bateman dual system and 
Pancharatnam's phase boils down to the ordinary Berry- 
Anandan phase. In Section V we take this observation 
over to the radial wave functions and then, by setting 
I = ±i, to the ID l.h.o.. Because the harmonic oscillator 
wave functions obtained in this way are constructed en- 
tirely from the fundamental system of solutions of Bate- 
man's dual system, the Berry-Anandan phase bears an 
imprint (or memory) of the original 2D system even after 
the reduction to the ID l.h.o. is performed. The geomet- 
ric phase thus obtained can be directly identified with 
the zero point energy of the ID l.h.o., and, in general, 
it is different from the usual Eq — Ml/ 2. This is in line 
with the results obtained in |]l7| . 

II. TIME-EVOLUTION AMPLITUDE (KERNEL) 
FOR BATEMAN'S DUAL SYSTEM 

In the following a decisive role will be played by the ma- 
trix elements of the time evolution operator U{th,ta) in 
the localized state basis 

(xb;4|Xa;ta) = {Kb\U{tb,ta)\:X.a) ■ (1) 

They are referred as time-evolution amplitudes, or simply 
kernels. Due to the fact that U{tb,ta) fulfills the time 
dependent Schrodinger equations 
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ih—Uitb,ta) ^HUitb,ta), 
Otb 

d 

lh—U{ta,tb) = -U{ta,tb)H ,tb>ta, (2) 

Otb 

the kernel satisfies the equations : 
d 

lh—{Xb;tb\Xa;ta) = H {-lhd:iCh,Xb) (xfe ; ifc |xa ; ia) , 

otb 
d 



ih—{xa;ta\xb;tb) 
otb 



(3) 



^ -TH'^ {-ihd^^,yib)T ^ {y.a;ta\yib\tb) ,tb>ta, 
with the initial condition 

lim (xa; ta|x6; tb) = 5(xq - Xb) . (4) 

Here T is the (anti-unitary) time reversal operator and 
W is the Hermitian-adjoint Hamiltonian (in most appli- 
cations H is both Hermitian and even under time reversal 
so f and T are usually omitted). An important observa- 
tion is that for quadratic Hamiltonians the kernel has a 
very simple form, namely 



(xb;tb|xQ;ta) = -F[ia,4]exp |^-S'c/[x]j , tb > ta ■ (5) 

The function F[ta,tb] is the so called fluctuation factor 
pof and is independent of both x^ and x^ ]l8| , p^ . The 
form (|^) is usually attributed to Feynman and Hibbs 
1^,^, but one may readily see that it is nothing but 
the kernel version of the celebrated WKB approximation 
(often referred to as the Van Vleck formula), which 
turns out to be an exact relation for quadratic Hamilto- 
nians. 

For a system with a time-independent Hamiltonian, the 
kernel reads 



(xb;4|xa;ta) = (xb|exp H{tb - ta) j |Xa) . (6) 

Inserting the resolution of unity (completeness relation) 

^ |'0m)(V'm| = 1 ■ 
m 

(here lipm) are orthonormal base kets at t = spanning 
the Hilbert space) into (^), we obtain that 

{Xb;tb\Xa;ta) = ^1pm{^b,tb)i^mi^a,ta) ■ (7) 

m 

Here we have identified ipm{^,t) = {x\ip,n{t))- The sym- 
bol * denotes usual complex conjugation. Note that 
■!/'m(x, t) and ?A^(x, f), obey first and second equation 
in (|3|), respectively. 

In the following we shall use the Feynman-Hibbs pre- 



scription 
system Jll 



for the the quantization of Bateman's dual 
L3[ . Our analysis will also reveal a host of sub- 



A. Lagrangian and classical equations of motion 

Bateman's dual model describes a 2D interacting system 
of damped-amplified harmonic oscillators. The corre- 
sponding Lagrangian reads [lffl-Q4j22] 



(8) 



L = raxy + {xy — xy) — nxy 

giving the (classical) equations of motion: 

mxci + -yxci + KXci = , 
mijci - jilci + Kyci ^ , 



(9) 



It is interesting to observe that the equation for x de- 
scribes the damped harmonic oscillator, while the equa- 
tion for y characterizes the amplified oscillator. In addi- 
tion, with appropriate initial conditions both systems are 
mutual mirror images. In this sense it may be sometimes 
helpful to think of y as describing an effective degree of 
freedom for the reservoir to which system with the x de- 
gree of freedom is coupled [ pT| , p"3[p^ . 

In the following it will be useful to work with the rotated 
variables fTJ]: xi ^ {x + y)/V2, X2 = {x - y)/V2. Then 



L = — (Xi - X2) + ^{XlX2 - X2X1) - -{X^ - X2) 



m , 

m 7 K 
-XX+-XAX--XX. 



(10) 



where we introduced the notation ab = gafi a°'b^ , aAb = 
e"^aabj3 and x" = {xi, X2) with the metric tensor gap — 
(o'3)q/3 (note also that e"^ = — £0/3)- The corresponding 
conjugate momenta read 

p = mx-i7crix. (11) 

In {xi,X2) coordinates the equations of motion read 

TO Xc/ + 7(71 Xc; + KXcl ^ (12) 

Notice that if u(<) is a solution of ( p^ so are (Tiu(t), 
cr3u(— t) and ia2u{—t). For the future reference it is use- 
ful to realize that the Wronskian is t independent (i.e. it 
is a time invariant of the system). Indeed, in our case 
the Wronskian has the form (to is arbitrary) : 



W{t) = W{ta)e^p 



dtTr 



to 



Vto / 



(13) 



Eq.(|13D is nothing but Liouville's theorem of a differential 
calculus applied to (|2h. 



B. Classical action 



tleties which are hidden in the seemingly clear relation Using the usual definition for the action: 
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S[x] ^ J dtL, 



we can write 
j-t, 

S'cz[x] = 

7 



dt 



m f d 



2 \dt 



l{xiXi — X2X2) — XiXi + X2X2 



-(X1X2 - X2X1) - -(Xi - X2] 

-^{xixi — X2X2)\tl — J at — (mx + 7(Tix + kx) 



Ccl{tb)±cl{tb) - Xc;(ia)xd(ia)] ■ 



(14) 



C. Fundamental system of solutions 

A fundamental system of solutions (i.e. a maximal sys- 
tem of linearly independent solutions) for Eq.(^2|) con- 
sists of four real 1x2 vectors {i = 1, 2, 3, 4). The rea- 
son why there are four independent solutions is, roughly 
speaking, a result of the fact that we have two boundary 
conditions for each index. Independence of solutions may 
be checked via the Wronskian, which has to be non-zero 
at least at one time t (actually the Wronskian is time 
independent here). In our case the Wronskian is the de- 
terminant of a 4 X 4 matrix: 



Wit) = Wito) 



Ui U2 U3 U4 
Ui U2 U3 U4 



(15) 



An important technical simplification may be achieved 
by realizing that we may always find such a fundamental 
system where two arbitrary solutions (say, U3 and U4) 
are set to zero at ta- This is due to the fact that in order 
to fulfill the boundary condition on x(<a) we need only 
two linearly independent vectors. Let us fix the following 
convention: U3 = vi and U4 = V2. Then the condition 
on the fundamental system may be rephrased as 

W{t) = W{ta) = \ui{ta)u2{ta)\ X |vi (t,)v2 (ta) | ^0- 

Note that if we had assumed the existence of a fundamen- 
tal system having three linearly independent solutions 
being zero at ta, the Wronskian would vanish identically. 
Any real solution of (^2|) might thus be written as 

Xcl{t) = QlUi(t) -t- a2U2{t) + (3lVi{t) + I32^2{t) , 

with ai and j3i being real numbers. Applying Cramer's 
rule, the solution Xc/(t) with two fixed points Xc/(ta) = Xq 
and Xci{tb) = Xf, reads 



Xc/(t) 



[Ui (t) Di + U2 [t) D2 + Vi [t) D3 + V2 (t) Di] 



UaVb 



(16) 



where Ua = \ui{ta)u2{ta)\, Vb = \-vi{tb)v2{tb)\ and 
Di 



Xa U2{ta) 
Xfc U2(4) Vi(tb) V2itb) 



Do 



Ui{ta) Xa 

ui(ifc) Xb vi(tb) V2(ih) 

Ul(ia) U2(<a) X^ 

Ui(4) U2(tb) Xb V2(4) 

Ul{ta) U2{ta) Xq 

Ui(4) U2{tb) Vi(tb) Xb 



Xa A U2(ia) X Vb , 
-Xa AUi{ta) X Vb, 



(17) 



An equivalent, and more useful, way of writing Xci{t) is 
to expand it in terms of Xq and Xb. After some algebra 
we get 



Xc;(<) 



where 



[xl Bi(t) + xl B2{t) + xl Bsit) + xl B4(t)] 



UaVb 



^'^'^-[Blit) 



Bl and are given by the determinant D, 



(18) 
(19) 



D = 



Ul{ta) U2(ta) 
U2{tb) U2(tb) Vi(tb) V2(4) 



= UaVb , (20) 



with i-th row substituted by {u\{t), U2{t), vl{t), vl{t)) or 
{ul{t),ul{t),vf{t),V2{t)) respectively. So for example: 



Ul{ta) U2{ta) 

ulit) uUt) vl{t) 

ul{tb) ul{tb) vf{tb) V^{tb) 



(21) 



As a result, the classical action iScJx] might be written 



as 

Sci fxl 



m 
2D 



xl XbBi(tb) + xl XbB2{tb) + xl XbB3(tb) 

-a;bXbB4(ib) - xl XaBi{ta) - xl XaB2(ia) 
xlxaB3ita)-xlxaB^{ta)\ . (22) 



An explicit representation of the action is 
m 



Sci [x] 



Bl{tb, 



2^ \ -{xlfB\{ta) + {xlfi3l{ta) 

+ {xlfBl{tb)-{xlfBl{tb) 

+xlxl (Blita) - Blita)) +xlxl (i?l(ib) 
+Xlxl {Blitb) - Blita)) - Xlxl (iilitb) + i3l{ta) 
+Xlxl (i3l{ta) - Blitb)) + xlxl (iilitb) + Blita) 

Using the basic properties of determinants it is possible 
to show now that both 5c/ [x] and Xc;(t) are independent 
of the choice of the fundamental system of solutions. We 
show this in Appendix B. 
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D. Fluctuation factor 

We can now take advantage of the Feynman-Hibbs obser- 
vation about the time-evolution ampHtude (ker- 
nel) for systems governed by quadratic Hamiltonians: 



{Kb;tb\Xa;ta) = F[ta,tb]cxp [ ^'^'dM 



tb > ta ■ 



(23) 



As remarked, the fluctuation factor F\ta,th] is indepen- 
dent of x„ and x^. In addition, from (23) follows that 



F[ta,tb] = (0;4|0;ta) = (0|C/(ifc, ia)|0) , 
and so for time independent H one has F[ta, tb] = F[tb — 

ta]. 

The most usual way of calculating the fluctuation factor 
is via the Van Vleck-Pauli-Morette determinant [E3-G7| : 



F[ta, tb 



det 



det 2 



2nn d4 ^ 



Zdeto 



2TTh did 



(24) 



The symbol det2(...) denotes 2x2 determinant. In 
( p^ ) we have also used the identities Pa — —^f^ and 
Pb — (one should take a little care when using the 
covariant and contravariant indices). Indices a and b are 
kept fixed throughout calculation. 

Actually ( p^ is correct only for sufficiently short elapsed 
times tb — ta as it was proved by Pauli [ p8| . In the 
general case, the determinant on the RHS of (|2^) will 
become infinite every time the classical (position space) 
orbit touches (or crosses) a caustic. A detailed examina- 
tion of quadratic systems reveals @,|l^ that (HI) remains 
valid even after passing through the caustic, provided we 
write the fluctuation factor as 



F[ta,tb] = 



det2 



i dpb 



(25) 



and insert a factor exp(— i7r/2) for every reduction of the 
rank of 1/ det2 (^d^ Sd/dx.^dx.'^^ at the caustic. Thus we 
have {tb > ta): 



of the classical path running from Xa to Xf, The form 
( p6| ) is due to Gutzwiller |3^ and the prescription ( [2^ ) 
is nothing but the connection formula for relating the 
kernels on both sides of the caustic in a continuous way 
To simplify the discussion we omit for a while the 
delicate issue of caustics assuming that the determinant 
in ( p^ ) is positive. We shall, however, return to it in 
Sections IV and V. 

Now we are ready to calculate the fluctuation factor. A 
little algebra gives us 



F[ta,tb] = [- [Bl{tb)Bl{tb) - Bl{ta)Bl{tb) 

+Bl{tb)Blita) - Bl{ta)Bl{ta) + Bl{ta)Bl{tb) 
+Bl{ta)Bl{ta) - Bl[tb)Bl[tb) - Bl{tb)Bl{ta) 




(27) 



Here we have used the equations of the motion and the 
fact that the Wronskian is time independent. Since the 
kernel is uniquely determined from the classical action, 
our argument on the uniqueness of the classical action 
implies that (x;,; tf,|xa; t^) does not depend on the choice 
of a fundamental system. Note also that due to the 
fact that F[ta,tb] = F[tb - ta], it follows from (|2^) that 

D{ta,tb)=D{tb-ta). 



III. WAVE FUNCTIONS FOR BATEMAN'S 
DUAL SYSTEM 

A. Wave functions il)n^i(r-,u,t) and xf}n.t{r,t) 

In order to calculate the wave function it is useful to 
rewrite the kernel in hyperbolic polar coordinates (r, u), 
with xi — rcoshu and X2 — rsinhu, and then apply the 
defining relation JlSt : 

{rb,Ub;tb\ra,Ua;ta) = ^ V'«,i(^fc: 4)V'i*; (^a, "a, ^a) , 



tb > ta- 



(28) 



Here we have used the symbol {*) instead of the usual 
complex conjugation symbol * - the need for this refine- 
ment will show up in the following. Invoking (|^), (|22' 
and (|2^) (see also Appendix C) we obtain: 



(X6;ib|x,;0 =e-' = "-''i^[ta,i6]exp(^-5,,[x]j . (26) 
Here na,b is the Morse (or Maslov) index |,0,|l9|,|9|-|3|] 



*The set of all points where the inverse of the Van Vleck- 
Pauli-Morette determinant vanishes is called a caustic. The 
Morse index then counts how many times the classical orbit 
crosses (or touches) the caustic when passing from the initial 
to the final position. In the literature, crossing points are 
often called focal or conjugate points. 
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{rb,ui,;ti,\ra,Ua;ta) 



m W 
2^V "D 



exp 



2Dh 



2rarb Bl{tb)cosh{Au) - 2rarh BKh) sinh(A 



m W 



'JH U2 , 21 B\{ta) 



B]{tb) WAN B^{tb) .WAN 

' cosh(Aw) sinh(Au) 



D 



D 



Au^ Ub- Ua, tb> ta, 

By observing that 

[Biitb)r-[Biitb)r = wD., 

we can put 



(29) 
(30) 



Blitb) 



D 



— cosn a , 



Bfih) fw . 
— ^-^ — - — \ — Sinn a . 
D \l D 

In Appendix D wc show that 

a{ta,tb) {ta ~ tb) + f3 . 



(31) 



Here T — ^ and /3 is a complex constant. It is also 
useful to denote the reduced oscillators frequency as 

— — . If not indicated otherwise, Q will be 

assumed to be real throughout. That is, we shall mostly 
be concerned with the under-damped case although oc- 
casionally a result can be taken over to the over-damped 
case. 

Eq.(^ll) allows to rewrite the kernel ( p^ in the following 
form 



{rb,Ub;tb\ra,Ua;ta) 

m fw 

■ exp 



2TTh \ D 

im IW 



im j dD 2 dD 2 
4Dn V f^^" ~ dh^'' 



X exp 



D 



raTb cosh(Ait — a) 



(32) 



This expression may be recast into a more suitable form 
if we apply the Laurent expansion pl|,p2[[i] 



(33) 



exp(zacosh(u)) = ^ (— 1)'/;(— 



/— — 00 



(/;(...) are the modified (or hyperbolic) Bessel functions) 
together with the addition theorem for the generalized 
Laguerre polynomials (Mehler's formula pi|,p^): 



11=0 



r(n + ; + i) 



(Z1Z26)" 



■ exp 



1-6 

For this purpose we set 



-b 



zi + 22 



1 - b 



h 2 



1-6 



(34) 



Zi = — yW 



P(taj n p(tb) 



with r'a; > 00 and 

p{t) = V{t)W 



dt Y 
W)J 



V{t) 



\ 



(u»(OAUj(i))' 



(35) 



t<3 



This allows to identify the parameter 6 appearing in (|3^ 
with b{tb) where b{t) reads 



bit) 



^t^/V{{) + ^p{t)-V{t) 

p(t) 



exp —i2 arcsin< 



Pit) 



where V{t) = u^it) A u^it) ^ vi(t) A V2(t). Note that 
p{ta) = Ua and that b{ta) = 1. 

The previous manipulations permit us to formulate the 
kernel in the desired form 

{rb,Ub]tb\ra,Ua]ta) ^ -} / 



TT — ^ T{n + 1 + 1) 



X [b*{ta)b{tb) 



\n+- 



w 



h y p{ta)p{tb) 



X exp 



2h 

m 
2h 



Pita) 

1 Pitb 

2 pitb) pitb) 

ipjta) VW 
2 pita) pita) 



Pitb) 



iranY 



^2 gi("a-tl6 + a(ta,ti,)) 



(36) 



^Because exp(i acosh(u)) is an analytic function of w - the 
only essential singularities are in ?i = ±00 - the Laurent ex- 
pansion (pSl) is well defined for any complex u. 



Because cosh(u) > sinh(u) then \xi\ > \x2\ and so we have 
automatically that > is a kinematic invariant. 



6 



Let us identify the wave function ipn^iir^Ujt) = 
(r, u\ipn^i{t)). Note that Eq.(|3^) immediately imphes that 

ij}^*\{r,u,t) cannot be associated with tp^ i{r,u,t). It is 
not difficult to see that this peculiar behavior goes into 
account of the seemingly harmless expansion (|3^). The 
point is that we have tacitly used the discrete (Laurent) 
expansion even if an alternative integral (continuous) ex- 
pansion was available This favoritism towards dis- 
crete Ps was deliberate (see also next Section). A care- 
ful analysis will reveal that the discreteness of I is not 
compatible with a unitary representation of the dynamic 
symmetry group of the theory. The remedy will be found 
in a self-adjoint extension of H and it will turn out that 
= tpn,iir,-u,-t). 

Now, from ( p^ ) and ( ^ ) we may deduce the wave func- 
tions 



,1 



ttV r(n + z + 1) VV 



l+l 



Pit) J 



i p{t) VW 



^2 g-i(„+rt-f ) 



V'i1(^,«,i) = \/e 



ttV T{n + l + l) VV 



X cxp 



2h 



i p{t) y/W 



^2 gi(«+rt- 



Obviously, neither ■!/)„^;(r, m, t) nor ^pl*](r,u,t) belong to 
ordinary Hilbert space because they cannot be normal- 
ized in the usual manner (they do not belong to the space 
of square integrable functions i'^). The latter observation 
is in agreement with Refs. ||ll|,|l^. and we shall com- 
ment more on this point in the next subsection. We note 
that the kernel (^) (and consequently the wave functions 
(^)) satisfies the time-dependent Schrodinger equation|^: 

■ifl-^ H{n, Ub)^ {n, Ub] tb\ra, Ua] ta) ^0, h > ta , 

where 



(37) 



H 



2m 



—pI + rn^n'^r'^ 



Lpu 



(38) 



1 

2m 



—fi^ I \- m il r 

Qj.2 J, Qj, j,2 Qy2 



ou 



The Hamiltonian ( |38| ) is the so called Bateman Hamilto- 
nia n W% - We now define the radial kernel {rb',tb\fa\ta)ni 
as 

The corresponding wave function ipn.iiT^t) = {r\ipn^i{t)) 
reads 



n 



m 



.14^1/4 



r(n + ; + i) VV f^pit) 

„2 1 



l+l 





2 p{t) pit) 



(39) 



It is simple to persuade oneself that ■)/'«,;(?', i) fulfils both 
the orthonormalization condition 

dr%lj'^*](r,t) 'ip„>^i(r,t) = S„„> , 
and the resolution of unity 

oc 
n=0 

Note that both the radial kernel and the wave function 
( ^ ) satisfy the time-dependent Schrodinger equation 

d - \ 

^^0^ Hlin) j {rb]tb\ra; ta)n,l = , h > ta , 



where 



Hi = 



2™ 



- 1 I 2o2 2 

- ] + m It r 



— iha[t) I — ihT I . 



(40) 



The term proportional to 1/r^ is analogous to the cen- 
trifugal barrier known from rotationally invariant sys- 
tems and so the quantum number I can be viewed as 
analog of the azimuthal quantum number. Note that, 
due to the structure of a{ta,tb), the term a(t) + F must 
be zero. 

Since the generalized Laguerre polynomials Ljj are de- 
fined for alH e C indicesQ, the wave functions ( p9| ) sat- 
isfy the time-dependent Schrodinger equation with the 



^Actually il)^*\{r,u,t) fulfills the time-reversed (time- 
dependent) Schrodinger equation, see Eq.(|4!l). 



**Analytic continuation for I with K(Z) = —1, —2, —3, ... is 
however required, see e.g. |33l 



7 



Hamiltonian ( pO| ) also for non-integer Vs. The key ob- 
servation then is that if we continue / to the values 
the Hamiltonian Hi describes the ID l.h.o.. If we make 
use of the rules connecting Hermite polynomials with the 
Z = ±i La guerre polynomials we may rewrite the 
continued radial wave functions in a simple form 



i'n^ir, t) 



1 



1 



22»+i Y n! r (n ■ 



hp{t) 



w- 



2n+l 



hp{t) 



r 



2 pit) pit) 



B 



1 



\/2hnnVl 
1 



[pi — imftxi] , 
[p2 ~ imflx2] , 



with 



(43) 



[A, A^ - [B, St] = 1, [A, B] = [A, St] = . (44) 
The Hamiltonian (BSf) can be then rewritten as 



H = hniA^'A - b'^b) + ihriA^B^ - ab) 
= 2hinc - TJ2) , (45) 

where we have made explicit the associated 50(2, 1) = 
SUil, 1) algebraic structure: 



^n,-^ir, t) 



1 



22" \ nl r I 



hpit) 



X [bit)r+-^H2n 



hpit) 



W ^ r 



X exp 



2h 





2 pit) pit) 



(41) 



In passing we mention that the quantum numbers n and I 
appearing in (|36|)-(^7|) have been seemingly independent. 
So far the only obvious restriction was that n > inte- 
gers. However, for a consistent probabilistic interpreta- 
tion (in the r variable) and analytical continuation (41), 
the wave function 'ipn,i{r,u,t) is required to be bounded 
for |r| < 00. Using the asymptotic expansion for Ll^{z) 
(see e.g. |l|,||): 

LU^)-^,^^.F,i^n,l + l;z) 



nlTH + 1) 

r(n + ; + i) 



1 



I 



i-z)" 
^ 0. 



where i_F'i(— n,/ -|- l;z) is the confluent hypergeometric 
series |2^. It is not difficult to see that the only allowed 
values of n and I at which ^n,iif,u,t) fulfills above re- 
quirements are those where: 

2n + Z + l = 0,±l,±2, ... , 1^0. (42) 

So 1^1 > 1 and n > 0. 

B. Meaning of quantum numbers n and I 

Let us now consider the meaning of the quantum numbers 
n and /. To do this we must first understand the algebraic 
structure of the Hamiltonian (^). In Refs. [ p]Jl3| the 
following ladder operators were introduced: 



-iA^A~B^B)\ 



J+ 

,h = \iA^A + B^B 



At St , J_ 



AB, 

+ 1): 



[J+,J-]=-2J3, [J3,J±] = ±J± 



(46) 



Here C is the only Casimir operator ( SUil, 1) has rank 
1). In addition, [C,H] = [J2,H] = 0. If one defines 



Ji = i(J+ + J_) 



J2 



(47) 



the more familiar SUil, 1) algebraic structure appears: 
with 



J? 



•^1 + 4 



(48) 



It is simple to check that 
1 



C 
J2 



1 



"2 , 2o2 2 

—Pu + m it r 



2h 



Pu ■ 



In this connection it is important to re cog nize that the 
system described by the Hamiltonian (^5|) is both con- 
servative and invariant under time reversal. Because the 
latter point has been treated in the literature in a some- 
how ambiguous fashion (cf. Ref . ||l^ ) , we discuss it in de- 
tail in Appendix E. We prove there that TCT^^ = C and 
T J2T-1 = J2 from which it follows that THT-^ = H. 
It is precisely this time-reversal issue which obscures the 
quantization of Bateman's system, bringing about many 
subtleties which are difficult to grasp without an explicit 
knowledge of the time-dependent wave functions (|37|). 



Now, the crucial observation is that although from (|47 
i.e. from the very definition, J2 appears to be Hermitian, 
(^) implies that it has a purely imaginary spectrum in 
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|'0n,i(i)) (and this holds for all t) . The root of this "patho- 
logical" behavior is in the non-existence of a unitary ir- 
reducible representation of SU (1,1) in which J2 would 
have at the same time a real and discrete spectrum 
Nevertheless, both the discreteness and complexness of 
J2 spectra are vital in our analysis since they bring dis- 
sipative features in the dy namics: this was also the case 
considered in Refs. 1 11 -1^], to which we are going to com- 
pare our results. Thus the usual unitary representations 
of SU{1,1) are clearly not useful for our purpose. On 
the other hand, by resorting to non-unitary representa- 
tions of SU(1, 1) (known as non-unitary principal series 
we lose the hermiticityP] of J2 and hence the 
spectral theorem along with the resolution of unity. 



Actually the situation is not so hopeless 
deed redefine the inner product |13 3; 



One may in- 
to get a unitary 
irreducible representation (known as complementary se- 
ries [ ^^ ) out of non-unitary principal series. This 
may be easily done when we notice, using (|2^), ( |37| ) and 
Schwinger's prescription (^), that the states i/'i*; are not 
simple complex conjugates of ipn,i because they fulfil the 
time-dependent Schrodinger equation 



(^ih ^ + Tijt (r, u)T-'^ (r, u,t) = 0, 



(49) 



with the (effectively) non-Hermitian Hamiltonian. Ac- 
cordingly, what we have loosely denoted in ( p8| ) as 
ipl^]{r,u,t) is actually 

{Tipnj{-t)\r,u) = {ipn.i{-t)\r, -u)* = ipn.i{r, -u, ~t) , 

(50) 

as can be also double-checked from the explicit form 
(p7|). For the sake of simplicity we use [T |-!/'„_;(t))]t — 

{Tijjn,i (t) I . Clearly, if J2 were Hermitian then ipl^] — ip* ; 
as one would expect. 

The above considerations have some important implica- 
tions. To see this, let us rewrite the kernel ( |36| ) by means 
of the states \^n,i{t))- 

= ^(''6,Mb|'0n,i(4))('?''0«j(-*a)ka,'"a) ■ (51) 



We can formally introduce the conjugation operation 
("bra vector") as (V'n,;(<)| = [T |^/i„,z(-i))]t. Then the 



t^To be precise, we should talk about self-adjointness rather 
than hermiticity, but we shall assume here and throughout 
that this ambiguity does not cause any harm in the present 
context. 



resolution of unity can be written in a deceptively simple 
form 



1 



(52) 



The price which has been paid for this simplicity is that 
we have endowed the Hilbert space with a new inner 
product. In this context two points should be stressed. 
First, under the new inner product |'!/'n,i(^)) has a finite 
(and positive) norm. Second, J2 is Hermitian with re- 
spect to this inner product. Indeed, integrating by parts 
we get 

(V'n.i(Ok2'0m,fc(O) 



drdu r 



(*) 



{J2 V'«,i(OIV'm,fc(i)> 



■>pm,kir,u,t) 



(53) 



In (^ we have applied ( |152D , (|50| ) together with the fact 
that the "surface" term is zero (if fc 7^ Z then integration 
w.r.t. the variable r gives zero ||2l[], ii k — I then the 
product ipn I ^m,i is u independent). 



Note that the (|5^) implies that \TJ2T — J2, or equiv- 
alently: T J^T-^ = J2. Indeed, 

(54) 

The spurious time irreversibility of J2 apparent in ( p^ is 
an obvious consequence of dealing with the non-unitary 
representation of SU{1, 1). This is understandable, since 
a mechanism which does not preserve the norm (dissi- 
pation) is inherently connected with time irreversibility. 
From a mathematical point of view, we can interpret the 
relation TJ2T~^ = J2 as a self-adjoint extension of J2 
in the space spanned by the \'4in,i{t)) vectors. It should 
be also clear that when J2 is time reversible (e.g., in the 
usual Hilbert space £^) it is also automatically Hermitian. 

A pivotal consequence of the above is that 



(55) 



Thus the time-evolution operator is unitary under the 
new inner product. It is this unitarity condition, intrin- 
sically built in the kernel formula ^ (and successively 
taken over by the Feynman-Hibbs prescription), which 
naturally leads to a "consistent" inner product intro- 
duced in a somehow intuitive manner in Refs. Jll| , p^ . 
From now on the modified inner product will be always 
tacitly assumed. 
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So far we have dealt with the pecuHar structure of the 
Hilbert space. To interpret the quantum numbers n, I 
labeUing the constituent states, we start with the obser- 
vation that from the expUcit form (37) one can readily 
construct the Hermitian operator C (commuting with J2) 
which is diagonalized by ip„^i{r,u,t). Indeed one may 
check that 



J2 4'n,l{r,U,t) 

C-ipn.iir,u,t) - 
Here 
C = C 



id, , I , , \ 
2^ V'«,i(^"7^) = ^^Vn,i[r,u,t) , 



2np 



(2n + / + 1) ^PnJir,u,t). 



(56) 



PL 

4p2 



i p f d 



c2C 



R{t)C R-^{t). 



(57) 



The unitary operator R{t) has the form 

R{t) = ^(^;^)exp(^CGA)exp(^CGi3) 



5(^;0 =exp {lir^ 



with 



and 



r ( ^ 



^ Ah p 



The reader may recognize in Ga and Gb the SU{1, 1) 
displacement operators (i.e., generalized coherent states 
generators) ||3^. Alternatively, one may view Ga and 
Gb as the single- mode squeeze operators One may 
also notice that J2 is nothing but the generator of two 
mode-squeeze |37 . Actually, the fact that there should 
be a close connection between SU {1,1) squeezed states 
and damped oscillators was firstly proposed in Ref. . 
Finally we should point out that in deriving (|5^) the 
relation pr = —ih [d/dr + l/2r) was used p^ . 

From Eq.(|5^) wc find that C is Hermitian whenever p 
is a real function and TC{t)T^^ — C{—t). The latter 
together with (|56|) implies that 



T|V„,i(0> = \^n+l,-li-t)) 



(58) 



and so in the static case (i.e. when p{t) — const, and 
V{t) = p sin^{nt), see also Section V) we have 



J2 ICi 



(59) 



Here (and throughout) the convention \ip^ ;) = l?/)^ ;(0)) 
is employed. Notice that in view of relations ( |152| ) and 
( ^ ) , the time reversed stationary states fulfill the condi- 
tion: 



n+l,-l/ 



(60) 



Using (p6D, { p7\) and (p9D, we can find the relation be- 
tween \'ipn,i{t)) and the stationary states The fol- 
lowing relation holds: 



(61) 



Thus the vectors \ipn,iit)) appearing in the spectral de- 
composition of the kernel ( p8[ ) have a tight connection 
with SU{1, 1) coherent states and, as we shall see soon, 
they describe indeed coherent states which (if expressed 
in the |r, u) representation) rotate in their position spread 
and/or pulsate in their width ("breathers"). 

However, before discussing other algebraic properties of 
the time-dependent states lipn.iit)), it is convenient to 
establish a connection with the results presented in Refs. 
[]ll| , P^ . Let us first denote by {|nyi, ns)} the set of eigen- 
statcs of A^A and B^fB. From (|^ follows that ua and 
ub are non-negative integers. Defining 



j ^ - {nA~ ub), 



-{ua + ub), (62) 



we can label the eigenstates of C and ( J3 — i) as |j, m) 
rather than |n^, ub)- As a result one may write 



C|j,m) = j |j,m) , J3|j,m)= TO-f- |j,m), (63) 



with two obvious conditions: |j| = 0, i, 1, |, . . . and m 



\]\Aj\ + \Aj\ + i 



The latter is compatible with wl 



Defining the vacuum state |0, 0) as the state fulfilling 
A |0, 0) — B |0, 0) = 0, then the discrete states \j, m) can 
be explicitly written as 



with 



|j, m)=c,-„ (At)^+™(i?t)™-^ |0,0) 



= [(j+m)!(m-j)r^/' 



(64) 



To relate the eigenstates of J2 with those of J3 above 
constructed, we may employ the following relation 

exp(6' Ji) J2 exp(— Ji) = — J2 cos6' — i J3 sin6' , 

with ^^ e C, to obtain that 

J2 = ±i exp Ji^ J3 exp (^=f^ Ji) ■ (65) 

The former implies that 
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,(±) 



I* 



(±)\ 



(66) 



Because C commutes both with J2 and Ji, we also have 

C|*S)=j'l*5^n)- (67) 

The relation coincides with the result found in . 
Similar identity was also realized in Ref. | |T^ . 



Comparing (pq) and (B7f) with (59) we can identify 



(+) 



2 T 2 ' 2 2 1 



(note, as TO > 0; / > 1) or equivalently 



Similar identification holds for 



(-) 



\rn,i) = * 

{I < —1) or equivalently 



(-) 



i^i — 1 — i / ' 
2 ' 2 ' 2 2 / 



I* 



(-)\ 



^1) 



(68) 



(69) 



(70) 



(71) 



Matching (|60|) with (|6|)-(|7l|), the states |*j;2> and 
l^tj- ^) can be related in a simple way, namely 

which can be interpreted as a continuation of 1^'^^') to 
negative to's. With ( |7^ ) we can check the consistency of 
the inner product defined above. Indeed, using (|66| ) we 
have 



(+) 



1 4- 1 1 

2 ^ 2 ' 2 



^ 2 ^ 2 ' 2 2 



— Snri '- 11' ■ 



(73) 



Now, to understand better the physical nature oi\ipni (t)) 
B 

left SU{1, 1) doublets, respectively. For instance, under 
an SU (1,1) rotation 



let us first note that ( | and ( I are right and 



Jv Jl 



B 



cosh| i sinh| \ f 



-i sinh-| cosh-| 



B 



M 



B 



(74) 



The transformation matrix M is then clearly an element 
of the SU{1, 1) group as M^g = gM^^ ; g = diag(l, -1). 



Analogous transformation rules hold also for "rotations" 
with respect to J2 and J3. 



Thus in terms of the ladder operators A and B, the states 



l*S> read 



^^^') = c 

^ l.ml ^3 



and so 



B^ ±A 

X exp[±7r/4(yl^B^ +^B)] |0,0) 

= {A^ + BT [B^ ^ Af\—' 
X exp[7r/4 (A^fit _^ ^5)] |o^ 0) , 



(75) 



with 



_ l)!2l'l-i 



-1/2 



It follows from (|68|) and ([70|) that ([75|) is true both for 
I > 1 and Z < — 1. Using the Baker-Campbell-Hausdorff 
relation one can find that 



exp[6' Jl] = exp 
X exp 



tan I - ) J+ 



tan I - I J_ 



exp 



2 Ig 1 cos - ) J3 



(76) 



(the formula ( |76| ) is an analog of the Gaussian decompo- 
sition well known from SO{3) group) and so ( [ts]) can be 
recast into a simple form 



IC,J = ^(^t + s)«(st + ^)i/i 



J (At)" (Bt)l'l-"-i 



-"-1 exp(J+) |0,0) 
(77) 



with 



A=^(A-Bt) B = ^(B-At) 
V2 V2 



1 



V2 



{A^ + B) , 



1 



V2 



(St + A) . 



(78) 



Since the canonical commutation relations are conserved 
by a similarity transformation ( |7^ ) , A and A^ are new an- 
nihilation and creation operators, respectively (the same 
holds true for B and B^), with a new vacuum state |0)) 
However, because the similarity transformation (^^ for 
rj = —iTr/2 is not a unitary transformation, A and A^ 
(and B, B^) are not Hermitian conjugates. This should 
not be surprising: we have already observed that the 
state {tpn i) does transform under a non-unitary represen- 
tation of SU{1, 1) and recognized this point as the origin 

of the "anomalous" behavior of ipn.i and '«/'„;■ However, 
it should be born in mind that under the inner product 
above introduced the A, A^ (and B, B^) are Hermitian 
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conjugates. Indeed, one may readily check that, for in- 
stance, ((0|A = ((AT 0|. 

In connection with Eq.([r7|) we should mention that |0)) = 
1/V2 [cxp( J+)] |0, 0) is a two-mode Glauber coherent 
state [0. The coherent state |0)) can be physically vi- 
sualized as a boson condensate or as a new vacuum asso- 
ciated with A and B operators. This suggests to us that 
the state li^n i) '^^^ ^le interpreted as an excited state of 
n type-A oscillators, and — n — 1) type-B oscillators 
with the coherent state type vacuum |0)) = |V'o - i)- 

Yet another, interesting interpretation of jV',^ ;) can be 
obtained if we employ the SU{1, 1) coherent states (i.e., 
generalized coherent states [Q). We demonstrate this in 
Appendix F. 



Using ( |77D it is now simple to interpret the quantum 
numbers n and / appearing in the time-dependent states 
\ipn,i{t))- Let us first observe that 



A^ (t) = R{t) At R-\t) ^aA^ -b B'f 
B\t) = R{t) Bt R-^{t) ^a*B^ - b* A^ 

('*' denotes a complex conjugation) where 



(79) 



As a matter of fact, the state 

\[CM = ^MKGA)eMKGB) |0)), 

is also a coherent state (it saturates the uncertainty re- 
lations) called two-mode squeezed state or (from 
rather historical reasons) (two-mode) two-photon coher- 
ent state 1 4^. Term squeeze (or squeezing) coined in Q 
reminds that although the dispersions of canonical vari- 
ables saturate the uncertainty relations their distribution 
over the phase space is distorted (or "squeezed" ) in such 
a way that the dispersion of one canonical variable is 
reduced at the cost of an increase in the dispersion of 
the canonically conjugated one. The concept of squeez- 
ing can be extended also to SU{1, 1) coherent states Q, 
however, the actual interpretation is in this case some- 
how less clear and states thus obtained do not seem to 
be particularly relevant in the present context. 

The physical meaning of |[CjCi^])) can be understood 
from the corresponding dispersions of Xi and pi\ 



2mn 

h 



{A-A^f)-{{A-A^)f 



2mn 



exp(-2C) , 



a= ( 1 + coshC--^- sinhC 



1-— - I sinhC+— - coshC. 



and |ap — \b\^ = 1. In a similar manner, we also find 

A{t) = R{t)AR-\t) ^a*A + b*B 
B{t) = R{t)BR-^{t) =aB + 6A. 



(80) 



Inverting Eqs.([79|) and (^0|) we immediately get the fol- 
lowing useful relations: 



a B 



V2A 

V2A^ = a* (a^ - B 



b* (a' 



blB^ 



a' 



^/2B = a* 
V2B^ = a (b^ - a) 



b ( B^ 



b* iaJ 



B 

FA 
A 

f B 



(81) 



By virtue of unitarity of R{t) the A{t), A (t) (and B(t), 

B\t)) are new annihilation and creation operators (albeit 
not Hermitian conjugates) with the vacuum state 

\[C,U]))^R(.t)m = s{C,t)\K,t])) . 

Hence, from (^^ we have 
IV'n.KO) =c„,;(A^t))"(B^i))"l— ' WCtm , (82) 



((Api)2) 



2 

hmfl 



{A + A^y)-{{A + A^)y 



exp(2C)+ (c/f^)'exp(-2C) 



(83) 



It is important to realize that ( p^j ) are obtained using 
( ^ ) together with the modified inner product. So, for 
example, the relation 

{{tA\c,tt]\ = {nA^\[C,^,-t]))]V ^^{{K,U]\B\ 

has to be employed. We may easily observe that if 
p — const, (i.e. S{^;t) — 1), the dispersions ( ^ sat- 
urate the Heisenberg uncertainty relations. As a result, 
the state |[C,i])) is indeed the squeezed coherent state 
with the squeeze parameter ^. If we evaluate the coherent 
state in the position representation we get a (minimum) 
wave packet specified by its half-width (stipulated via the 
Xi dispersion) and by the mean position (stipulated via 
mean of Xi) |^2|. If p 7^ const, (i.e. S{£,;t) ^ 1) we do not 
have any more a minimum uncertainty packet - uncer- 
tainty product does not stay fi,^/4 any more. Apart from 
the time dependence gained through C, (which would still 
allow for the minimum uncertainty wave packet) there is 
an additional contribution in the dispersion of pi as it 
can be directly observed from (|8^) . So depending on the 
time behavior of C, the wave packet width now oscillates 
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and/or spreads and this, in turn, "deflects" the disper- 
sion of pi from that of minimum uncertainty product. 

It follows that Eq. ( |8^ ) allows for a simple physical expli- 
cation of (r, u|^„^;(i)): the state (r, u|-i/;„_;(i)) describes 
the excited state of n type-A oscillators, and (|Z| — n — 1) 
type-B oscillators with the vacuum state (r, u| [C, ^, i])). 
The vacuum state itself is represented by the static wave 
packet which pulsates (and/or spreads) in its width. 
The excitations then may be understood as 2D Galilean 
boosts combined with SU{1, 1) rotations. 



IV. GEOMETRIC PHASE FOR BATEMAN'S 
DUAL SYSTEM 

A. Intermezzo: geometric phase 

Since in the following an important role will be played by 
geometric phases, we give here a very brief introduction 
to this subject. 

One has first to recall that a Hilbert space 7i is a line 
bundle over the projective space V, i.e. the equivalence 
class of all vectors that differ by a multiplication with a 
complex number. We shall denote a generic element of V 
as IV')- The inner product on 7i naturally endows V with 
two important geometric structures: a metric p3| , 3 

ds^^\\dm\''~mdm\ m 

and a C/(l) connection (Berry connection H,E5l) 



A = i{^\d\^)^ 



(85) 



When a point evolves on V along a closed loop, say 7, 
the total phase change (ptot of IV") on Ti, consist of two 
contributions: the dynamical part 



f^dyn 



-h- 



dt{m\H\m). 



(86) 



(with T being the time period at which the system tra- 
verses the whole loop 7) , and the geometric part (Berry- 
Anandan phase|^ 



= (V'(0)|V(T))exp(^z dtm)\i-m)) 

dtm)Vj^m)) 



exp I 







exp j A 



(87) 



So the Berry-Anandan phase ipBA may be geometrically 
understood as (an)holonomy with respect to the natural 
(Berry's) connection on the projective space V |Q. It is 
needless to say that all above relations are meant to be 
valid only for states |^/'(i)) being normalized to unity. In 
general, corresponding division of the state norm must 
be invoked. As our states 4'n,i{r-,u,t) are not normal- 
ized, this should be kept in mind in the following where, 
for clarity of notation, we shall omit this normalization 
factor. 

Actually, the geometric phase can also be defined for open 
paths. In this case the phase is usually referred to as Pan- 
charatnam's phase [Q. The trick is that any path on V 
can be closed by joining cndpoints with a geodesic con- 
structed with respect to the metric (0). The geometric 
phase of the loop thus constructed is then defined to be 
equal to the geometric phase associated to the open path. 
So, if 7o is an open path on V and 7^ is the corresponding 
geodesic on V, then the associated Pancharatnam phase 
(f)p reads 




d , ~ 



exp(^iy^ dt {^it)\i-\t^{t)) 
{i^{U)\^{tf)) exp (^i ^ ' dt m)\ij^Wt)) 



Here we have used the fact that parallel transport along 
a geodesic does not bring any anholonomy. Note that cjjp 
is well defined only if the endpoints are not orthogonal. 
It should be also clear that 0p is defined only modulo 
27r. 



"It would be perhaps more correct to call the geometric 
phase presented in this paper as Anandan's phase or 
Aharonov-Anandan's pha se [[43l . However, due to a historical 
reasons, the Berry phase [)46|| is usually taken as a synonym 
for any cyclic geometric phase. We feel that the name Berry- 
Anandan phase is a suitable compromise between rigor and 
tradition. 



B. Exact geometric phase in Bateman's system 

To find the geometric phase for Bateman's system we 
firstly compute the dynamical phase 

/■*/ d 

(t>dyn = - / dt {lpn,l{t)\i-^\lpn,l{t)) 

= -;i-iy^ dt {^/M,i{t)\H\^nAt)) ■ 

Using (13), (ll) and (|^) we find that 
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+ in {tf~ti) 



(89) 



In order to proceed further it is convenient to define tlie 
complex number 



z{t) = 



Pit) 



'1 - 



m 
Pit) 



Notice tliat zz = I. When evolving in the time interval 
(ti,tf), z{t) traverses a curve 7 in the Gaussian plane. 
The index of the curve 7 (i.e. the number of revolutions 
around the origin) is then defined as Es]: 



ind7 = 



1 



dz 
z 



(90) 



Consequently the dynamic phase can be rewritten in the 
form 



(t>dyn - -(2n + Z+ 1) 



dt 



W ' 2Vw 



{2n + l + l)(7rind7) + iTl{tf - t,) . 



n) 



Using (|90|) we may write also the total phase in a fairly 
compact form, indeed 

(t>tot = arg {{ipn,i{U)\ipn.iitf))} 

1) ' ~ 



= — (2ri + ^ + 1) ^arcsm^ 
- 7:ntj+ iTl{tf-U) 



Pitf) 



p{U) 



= -(2n + I + l)(27rind7) - - mj + iTl{tf - U) , 

(92) 

With riij being the Morse index of the classical trajec- 
tory running between x.^ and xj. 

At this stage a remark should be added. The fact that we 
get an imaginary piece both in (ptot and (t>dyn should not 
be surprising because we work with the modified inner 
product. Notice that the "troublesome" contribution in 
(ptot and (pdyn correctly flips sign if one passes from ■ .) 

to 



Substituting (|9l|) and (|9^) into the equation (| 
Pancharatnam's phase 



we get 



bp = I2n + I + 1) [ ^ dt( — ^ 

4 \8pVw 

- {2n + l + l) (7rind7) 

TT 



2VW 



(93) 



Note that because p, W and V are solely constructed out 
of solutions of the classical equations of motion, 0p is 
manifestly h independent. 



The meaning of the terms on the RHS of ( |9^ ) can be 
easily understood. Let us assemble the first two pieces in 
Eq.(B3h together. Using the fact that 



P 



n 



exp(-2C) , 



(see explanation below Eq.(p7D) together with ( |35| ) we 
obtain 

(j)p ^ {2n + l + l)- J dt (exp(-2C) + exp(2C) 



+ exp(-2C) {C/n)j+-n,j 



{2n + l + l) I dt' ^ ' ' ' ' ' 



hm 



TT 



(94) 



The index "i" is either 1 or 2 (it really does not mat- 
ter as dispersions are symmetric, one may write also the 
symmetrized version with the prefactor i, however, some 
care should be taken as we do not have Euclidean scalar 
product and so, for instance, (^2)^ = —P2P^)- It is impor- 
tant to recognize that (...) in ( p4[ ) represents the mean 
value with respect to the ground state. 

Thus (pp is a collection of three contributions; overall 
ground-state fluctuations of p and x gathered during the 
time period tf — ti and the Morse index. While the first 
two are basic characteristics of the ground-state wave 
packet, the Morse index contribution, on the other hand, 
reflects the geometrical features of the path traversed by 
the ground-state wave packet in the configuration-space. 
As explained before, its presence is inevitable for provid- 
ing a correct analytical continuation of the Feynman- 
Hibbs kernel prescription around the focal points. 

The above considerations show that the properties oi (pp 
are basically encoded in the structure and in the time 
dependence of the ground state. This intertwining of the 
ground state with geometric phase will be of a crucial 
importance in the following. 

An interesting question which one can raise in the present 
context is how the non-abelian (Wilczek-Zee) geometric 
phase looks in Bateman's system and to what extent 
it influences the presented results. This is definitely a 
challenging task as there does not exist at present any for- 
mulation of non-abelian geometric phases outside of the 
scope of adiabatic approximation (i.e. geometr ic p hases 
pioneered by Berry). As the geometric phases ([87|), ( |8^ ) 
and (^3|) are rather Aharonov-Anandan type, such an 
extension would be of a particular interest. We intend to 
investigate this question in the future work. 
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V. THE GROUND STATE OF THE ID LINEAR 
HARMONIC OSCILLATOR 

A. Berry— Anandan phase 

An interesting implication of ( p3| ) arises when we turn 
our attention to the special case of / — — i. From Sec- 
tion III we know that such a choice corresponds to the 
ID l.h.o.. In addition, from <^ and (|l|) we see that 
the following relation between the ID l.h.o. wave func- 
tion il?^°(r^ t) = ip„ _i (r, t) and Bateman's system wave 
function tljnd{r,u,t) holds: 



(95) 



The important point about ° (r, t) is that it is con- 
structed exclusively from the fundamental system of solu- 
tions corresponding to Bateman's dual system, and not, 
as one could expect, from fundamental system of solu- 
tions of the ID l.h.o.. In addition, it should be noted 
that the value u = —Ft + 13/2 entering (p5| ) is nothing 
but the solution of the classical equations of the mo- 
tion: u = dH/dpu, Pu = —dH/du with the classical 
Pu = J2 = and Bateman's dual Hamiltonian H. One 
may naturally wonder then whether some measurable in- 
formation about the original system could be tracked 
down in ■(/'^''"(r, t). We shall see that this is indeed the 
case: a memory of the underlying Bateman's system is 
imprinted in the ground state energy of the reduced sys- 
tem. 

Let us look at the geometric phase of the ID l.h.o. "in- 
herited" through the reduction (|9^) . Being r the mutual 
period of both p and V, then the following chain of rea- 
sonings holds: 

'^n°{r,'r) = \/^V'„- !(?-,", T)U=;3/2-Tr 
= (e'^-* V„,-i(r,u,0)) \u=p/2-rr 

= ,4^-H/o(.r(.)|.^_.|.r(*))^^,.^^^^^_ (96) 

On the other hand, because of the periodicity of p and V, 
the wave function ipl!/°{r,t) must be r-periodic as well. 
This implies that V'i''" (^, 0) = e'^''"^ 'ilj'^°{r,T), where m 
is an arbitrary integer. Using the fact that (pBA is defined 
modulo 2'K we can write 



(i)l^-i lVC°(i)) dt = h (2^n - <Pba) ■ (97) 



In particular, if ipl^°{r,t) are eigenstates of H_i, then 
we get the quantized energy spectrum: 



In the usual semiclassical treatment the presence of the 
Berry- Anandan phase modifies the energy spectrum via 
the Bohr-Sommerfeld quantization condition [^,0 . In 
the case of the simple ID l.h.o. (not the one obtained 
from the Bateman's system after reduction), the Berry- 
Anandan phase materializes only due to the Morse index 
contribution, i.e. (fiBA = —t^ na.b/2. When {tt — ta) — 
T — 27r/f7, the Morse index is simple [|l9|: Ua^b = 2 and 
one recovers the standard relation: 



1 



(99) 



In this respect the result (98) might seem rather pecu- 
liar, especially in the case when 4>ba The fact that 
4'BA can indeed be different from — tt will be explicitly il- 
lustrated in the following subsection. However, the basic 
reason for this to happen is not difficult to understand: 
the Morse index of the underlying Bateman's system (at 
u = —Tt ) is not necessarily equal to the Morse index of 
the corresponding ID l.h.o. at the same elapsed time r. 

Let us add two more comments at the end. Firstly, the 
foregoing analysis can be naturally extended on the case 
I = + i, but instead of doing that we may directly observe 
from ffl) that 



(100) 



so the ID l.h.o. states obtained from 1 (r, t) describe 
the higher energy states than ip'^°{r,t), and thus, for 
instance, the ground state comes entirely from the I — —-^ 
case. Secondly, it may happen that p and V can have 
more than one common period. Such periods can give rise 
to (generally) different (not mod(27r)) Berry-Anandan 
phases. This situation is fairly standard in many systems 
(see e.g., [Q) and the corresponding phases have as a 
rule different physical consequences. 



B. Practical example - stationary states 

To elucidate the previous analysis, we consider here an 
explicit example in which the following fundamental sys- 
tem of solutions is chosen: 



e: 



Iho 



{2nn - (I)ab) ■ 



(98) 



ulit) 

ulit) 
ulit) 
vl{t) 

vlit) 
vlit) 



V2co8{nt) cosh (Ft) , 
-\/2cos{nt) sinh(Ft) , 
V2cos{nt) sinh(rt), 
-V2cos{nt) cosh(Fi), 
V2sm{nt) cosh(Fi), 
-V2sin(m) sinh(Ft) , 
V2sm{nt) sinh(Ft) , 
-\/2sin(m) cosh (Ft) . 



(101) 
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The Wronskian is 



The explicit form of the wave function is then 



W = 



^2 

-V2 

\/2r V20 

-^/2T 



(102) 



and the determinant D 

D ^ Asin^ n{h - ta) . 
As a result we have 

B\{t) = 4 cosh[r(i - ta)] sin[f7(4 - t)] s\n[n{U - ^a)] , 
Bl{t) = 4 cosh[r(if, - t)] sm[Q.{t - ta)] sm[n{tb - to)] , 
Bl{t) = 4 sinh[r(tQ - t)] sin[f^(if, - t)] sin[f7(4 - ta)] . 

The classical action has the form 



Sri — 



2 sin [n{tb-ta)] 

'^fa^b cosh [Ufj — U 



{{rl+rl)cos\^l%^ta)] 

r(ifc-t,)]}, (103) 



and we find that the fluctuation factor reads 



F[taM 



n 



2Trh I sin [n{tb — ta) 



(104) 



Eqs.(103) and (104) lead to the kernel: 



{rb,Ub;tb\ra,Ua;ta) 



n 



2Trh I sin [n{tb — ta) 



X exp 



2hsm {n{tb -ta)) 
2rarb cosh (Am - r{tb - ta))} 



{{rl+rl) cos {n{tb~ta)) 



(105) 



Note that the kernel is indeed independent of the funda- 
mental system of solutions. One may check that the ker- 
nel ( |105D satisfies the time-dependent Schrodinger equa- 
tion IM). 



We now rewrite the kernel applying the expansion (36). 
Remembering that there is an absolute value in (105) and 
employing the fact that 

V{t) = 2sin2(m), 
P{t) = 2 , 



b{t) 



we get 

{rb,Ub;tb\ra,Ua;ta) = - 



frnny+' 



T{n + l + l) V ft / 



X Lt 



^ g-if2(2n+i+l)(t6-t„) g-iAu-ir(t6-t„) 



(106) 



il)n,i{r,u,t) 



TrT{n + l + l) \ fi I 
The radial wave function is on the other hand: 



V'n,i('", t) 



I n\ I m\/?2 

7rr(n + ; + 1) I h 



f ^ 2 g 2h ' 



(108) 



Note that ipn,i{'r,t) is an eigenstatc of Hi. In passing it 
is also interesting to consider the reduced wave function 
\'4>)Hj which is obtained from the full one through the 
formula = exp[— ri9„] This amounts to sub- 

stitute u with (u — Tt) into the total wave function (107), 
which then becomes 



■ip„,i{r,u-Tt,t) 



X r e 2" Li 



7 I 



m 



1+^ 



7rT{n + l + l) 



j,2 \ ^-lu g-jO(2n+i+l)t 



(109) 



It is easy to verify that it satisfies the reduced 
Schrodinger equation: 

i9tVn,i('^,"-rt,i) - 2r!C^/'„j(r,u-n,t), (110) 

This result will be particularly important in the following 
paper of this series. 

Because V is periodic with fundamental period r = 7r/fi 
and because V'n.-i is an eigenstate of H_i, the energy 
spectrum of the related ID l.h.o. is done by the prescrip- 
tion (|9^). The corresponding ground-state energy can 
be calculated from (p3|) and (98). We obtain 



Iho 



(pBA "-Q,6 

-h = rvK 



(111) 

2r ^ ' 



We see therefore that the fundamental system (101) re- 
fiects the dynamics of the underlying Bateman's dual sys- 
tem in E\^° only via the Morse index Ua^. To find na,b. 



we first define a Lagrangian manifold [glQ]: L ~ {x, p = 
dSci/dx}, where Sd is the action taken as a function of 
the end point x {— x^). For quadratic actions the La- 
grangian manifold is clearly n-dimensional plane in the 
2rt-dimensional phase space. Because the starting point 
of our analysis was the (configuration-space) kernel we 
are primarily interested in orbits for which the initial and 
final positions are given. Of course, if the initial position 
and momenta were given, then we would have a unique 
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orbit, but since instead we have initial and final posi- 
tions, it is not clear that any such orbit exists, and if 
does, whether it is unique. Let us therefore consider the 
time evolution of L. As time progresses, the Lagrangian 
manifold evolves (foliates) the phase space^. If we de- 
note a point of the initial-time Lagrangian manifold as 
Za = (Pa,Xa) (notc that set of all such point forms the 
Cauchy data for the Hamiltonian dynamics) then due to 
quadratic nature of H the final-time point Zf, is related 
with Za via linear canonical transformation (symplectic 
matrix): = Sz^ , where § = §(t£,,tf,). Taking only the 
X part of zt, we may write 



Xfe. 



(112) 



Clearly, if §i(ta,ib) were invertible (i.e., if det(§i) 7^ 0) 
then for given points Xf, and Xa would exist only one 
Pa and consequently only one classical orbit would run 
between x;, and x;,. If however det(§i) = 0, then ei- 
ther none or infinitely many solutions may be obtained, 
depending on the ranks of Si and the corresponding aug- 
ment ma trix. So if det(Si) = and Xa and x^ are such 
that Eq.(112) is not satisfied, then there are no orbits 
arriving at Xf,. If, however, Xa and x^ satisfy Eg. ( 112 ), 
then there is an n-dimensional infinity of initial momenta 
which maps onto x^ and thus density of particles is infi- 
nite at Xf, - orbits are focused in the configuration space 
(the situation is schematically depicted in Fig.j|]). 




FIG. 1. Evolution of a Lagrangian manifold schematically 
depicted at three distinct times. Note that due to the initial 
condition (4) we have La = {xa,Pa = anything}, note also 
that the Lagrangian manifold Lj, is responsible for a caustic at 
Xfc. Five possible orbits are shown. 

It is not difficult to see that the precarious points x;, 
in which dct(Si) = are precisely the focal (conjugate) 



^^Due to its very definition, Lagrangian manifolds transform 
under canonical transformations into other Lagrangian man- 
ifolds Q. So namely Lagrangian manifolds evolve into other 
Lagrangian manifolds under (Hamiltonian) time evolution. 



points. Indeed, 



= det(§i) = deta 



dcto 



ax°x 



When passing through a caustic, det(§i) may change 
the sign depending on the rank of det(§i (ia, ^fc)) at the 
caustic. The caustic is said to have multiplicity k if the 
rank of det(Si(tQ, i;,)) is k. The change of sign will di- 
rectly influence the form of the fiuctuation factor F[ta, tb] 
(which is basically the square root of 1/ det[Si (ta, it)]) 
as the correct branch cut must be chosen. The phase of 
F[ta, tfc] can be consistently prescribed demanding conti- 
nuity in the kernel js). It turns out that a phase factor 
exp(— ifc7r/2) must appear when passing a caustic of mul- 
tiplicity k. As more caustics are passed along an orbit, 
the phase factor will accumulate. The Morse index Ua^b 
appearing in (26) then simply counts caustics (including 
their multiplicity) encountered by an orbit passing from 
the initial to the final Lagrangian manifold. 

Let us analyze our particular situation. Using the Hamil- 
ton equations of motion we easily get the following solu- 
tion for z = (pi,P2tX^ ,x^) (to be specific we work here 
in {x^ ,pi,p2) phase-space coordinates) 



= exp 
— cxp 
X exp 



and so 



{t - to) 

(t-ta) 
{t-ta) 



—mfl'^as 







-T^-CTi 
2m ^ 



sin[r^(tb - ^g)] 



CT3 . 



(113) 



(114) 



To g et the Morse index for the Berry-Anandan phase 
( |1 1 iD , we simply notice that during the time interval 
(0, r) orbits pass one caustic at the conjugate time r = 
n/n, with multiplicity 2 (see Figs.||,|,|). The Morse 
index is Ua.b = 2 and the ground state of the ID l.h.o. 
Eq'^° — hfl. Of course, if r is a period for V, 2t is also 
a period and any integer (positive or negative) times t 
is likewise. It is clear then that the Morse index based 
on any such period should be the same because V (and 
hence the wave function) cannot distinguish whether the 
fundamental period alone or its integer multiples are in 
use. Indeed, for instance, for the time interval (0, 2t) or- 
bits pass two caustics, at the conjugate times r = 'l/^ 
and r = 27r/0, both with multiplicity 2 (see Figs.||,|,||). 
The Morse index is in such a case = 4 and we get 
again EH"" = Kl. 



We see then that although our prescription does not al- 
low to pinpoint i?""' precisely (the ground states are de- 
fined mod(Sr2)), there is no way how to bring the usual 
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I2 



FIG. 2. Evolution of a Lagrangian manifold in {xi,pi) 
phase space. Three orbits with different values of pa are plot- 
ted. The Lagrangian manifold is shown at five different times. 
Note that times t — n/fl and t — 2-K/fl correspond to the con- 
jugate times. This and the following plots are done for k = 40, 
7 = 1.2, m = l. 




FIG. 3. The same plot as Fig. 2 in the {xi,p2) phase space. 



fraction 1/2 into the result. The factor 2 is a "mem- 
ory" of the underlying 2D system. The fact that this 
may happen is not difhcult to understand. The usual 
derivations of the ground state energy hinge either on the 
Heisenberg-Weyl algebra or directly on the Schrodinger 
equation and thus have only local character while Berry- 
Anandan phase is (non-local) global characteristics of a 
system. We may thus conclude that the ground state 
of the reduced system can generally change according to 
the global properties of the original underlying system on 
which the reduction is performed. 

Actually the above analysis is still not the whole story. 
In the paper to follow we will show that there is yet an- 
other - not so far considered - contribution to the ground 
state energy, reflecting the dissipative nature of the sys- 
tem when working with the SU{1, 1) non-unitary repre- 
sentation. Such a contribution will manifest itself in the 



FIG. 4. The corresponding {xi,X2) configuration-space or- 
bits. Focal points occur in Xa and x.b at times t = n/Q and 
t — 2Tr/Q, respectively. Irrespective of the initial-time mo- 
menta, all orbits reach the focal points at the same time. 



form of an additional phase factor - "dissipative" phase 
(see also 0i)- 



VI. CONCLUSIONS 

In this paper we have studied the quantization of Bate- 
man's dual system of damped-antidamped harmonic os- 
cillators [ p^pT| by using the Feynman-Hibbs kernel for- 
mula. It has been known for some time that Bateman's 
system is difficult to quantize due to a multitude of con- 
ceptual problems (ll-13|. In order to address some of 
these problems and to improve our intuition for the com- 
plications involved, we have found it very convenient that 
classical mechanics lies manifestly at the heart of the 
Feynman-Hibbs prescription. Thanks to the fact that 
Bateman's dual Hamiltonian (Lagrangian) is quadratic, 
the kernel is fully expressible in terms of the fundamen- 
tal system of solutions of classical equations of motion 
and in addition, it is independent of the choice of such 
classical solutions. 



Using the spectral decomposition of the time-evolution 
amplitude we have been able to calculate the full wave- 
functions, fulfilling the time-dependent Scrodingcr equa- 
tion. This helped us to understand the reported con- 
troversy in the quantization of Bateman's system. We 
have shown that the above inconsistency has its origin 
in two interrelated issues: apparent non-hermiticity of 
the Hamiltonian and oddness of J2 under time reversal 
(and thus time irreversibility of the Hamiltonian). We 
have argued that both these "pathologies" are the con- 
sequence of a single fact, namely that what has been in- 
variably used in the literature (mostly implicitly) was the 
non-unitary irreducible representation of the SU{1, 1) 
dynamic group. The fact is that from the representa- 
tion theory it follows that there is no SU{1, 1) unitary 
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irreducible representation in which J2 would have at the 
same time a real and discrete spectrum p3|] . There are 
then two possible solutions. We may work with the uni- 
tary representation of SU{1, 1) ~ Bateman's dual system 
will then be free of pathologies but not particularly in- 
teresting from the physical point of view. The interest- 
ing features, e.g., dissipation along with time irreversibil- 
ity, enter the scene precisely when the non-unitary ir- 
reducible representation of SU{1, 1) is used. To treat 
the latter situation mathematically we have to face the 
non-hermiticity of J2 ■ The remedy naturally arises from 
the Feynman-Hibbs prescription and is based on a re- 
definition of the inner product: we have illustrated the 
mathematical and logical consistency of such a proce- 
dure. The relation with existing results ||Tl|-[r^ was es- 
tablished and the corresponding underlying SU{1, 1) co- 
herent state structure of quantum states was discussed 
in detail. The reader may contrast our approach with a 
somewhat more customary treatment of non-Hermitian 
quantum systems by means of resonant or Gamow states 

Although the kernel is invariant under the choice of the 
fundamental system of classical solutions, this is not 
the case for the wave functions. For one-dimensional 



quadratic systems it has been argued l20| that various 
choices of fundamental solutions correspond merely to 
different unitary transformations. However, this shows 
up not to be correct in our case. In fact, the basic fea- 
ture of Bateman's dual system is that states are unitary 
inequivalent under SU{1,1) symmetry. As a result, the 
wave functions and thus the geometric (Pancharatnam) 
phases are found to depend on the choice of the funda- 
mental system of classical solutions in a non-trivial way. 
It is also worthwhile stressing that the geometric phases 
here obtained are given only in terms of the parameters 
of classical solutions and are thus manifestly h indepen- 
dent (as should be expected from phases with entirely 
geometric origin). 

The crucial observation made in the course of our analy- 
sis is that when we analytically continue the "azimuthal" 
quantum number / to ±1/2, the reduced (or radial) wave 
functions formally fulfil the time-dependent Schrodinger 
equation for the one-dimensional harmonic oscillator. 
We have found that the geometric phase of the ID l.h.o. 
obtained via such a reduction is not equal to the ex- 
pected Berry phase , but bears a memory of the clas- 
sical motion of the original Bateman dual system. This 
"shadow" of the underlying 2D system originates from 
three sources: overall ground-state dispersions of p and 
X gathered during the period of evolution and the Morse 
index contribution. It should be noted that the built-in 
memory in the geometric phase is not a new idea, having 
been used in various contexts as the Born-Oppenheimer 
approximation [5^ and the dynamic quantum Zcno effect 



A remarkable feature of the reduction procedure is that 
it allows us to find a relation between the ground-state 
energy of the ID l.h.o. obtained after reduction and the 
Pancharatnam phase. To put some flesh on the bones 
and to demonstrate the mechanism at hand, we have re- 
sorted to a specific fundamental system of solutions. The 
corresponding wave functions then proved to be, after re- 
duction to those of a simple ID l.h.o., energy eigenstates 
and periodic with the period t, the latter being connected 
with the reduced frequency of the original system. The 
reduced geometric phase - Berry- Anandan phase - was 
then found to be directly proportional to the ground- 
state energy of the ID l.h.o.. It was shown that the 
ground-state energy is controlled by the Morse index af- 
filiated with Bateman's dual system (not with the ID 
l.h.o. itself !). Finding a Lagrangian manifold and fol- 
lowing its evolution in phase-space, we were able to track 
down the number of focal points in the interval (0, t) and 
hence to identify the Morse index. It turned out that the 
ground-state energy thus acquired is different from the 
usual Eq = hfl/2. 

Finally, we remark that the reader may find some res- 
onance of the method presented here with the results 
p^-p6| which suggest that the quantum mechanical en- 
ergy spectrum can be determined from purely classical 
quantities such as lengths and stability indices of the pe- 
riodic orbits alone. Whether or not this formal similarity 
can go any further is definitely a challenging question 
which is being investigated by the authors |l7|]. 
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APPENDIX A 

It is useful to list some of the expressions of Section II as 
they look in {x,y) coordinates. The Lagrangian reads: 

7, 



L = mxy + — {xy — xy) — nxy 

m , . 7 . K 

= — XX H X A X XX , 

2 2 2 ' 



(115) 



where x°' — {x,y) with the metric tensor gap = {cFi)af3- 
The canonical momenta are 



1 

P = TOX - -7(T3X . 



(116) 



and the classical equations of motion can be written as 

TO XcJ + 70-3 Xc; + K Xc/ = . (117) 
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Notice that if u(t) is a solution of (117) then a^u^t), 
(Tiu(— i) and ia2U.{—t) are also solutions. 
The Wronskian for the (x, y) system is 



W{t) = W{to) 



and the action reads 



exp(-£d.Tr(^a3)), 



(118) 



Sci [x] 



dt 



dt 



{xy + xy) -xy-xy 



+ ^{xy~ xy) - ^xy - '^xy 



m 
m 



^{xy + xy)\l'' 



dt — (rnx + 7(T3X + kx) 



[^cl{h)^cl{tb) - Xc;(ta)Xc;(ta)] ■ 



(119) 



and more explicitly 



m 



Sci\x\ = — 



2^ L-(^^)'^i(^")-(^')'^2(*'^) 

+ {xlf Blih) + {xlf Blih) 

-Xlxl (B\{ta) + Blita)) +Xlxl (i?^ (t,) + ijl (t,) 

+xlxl (bHu) - i?|(t,)) + x\xl (b\{U) - B|(t,) 
Finally, the fluctuation factor in (x, ?/) coordinates is 



m 






-Bl{ta)Bl{h) 


-B?(tb)i?i(t,)^ 


V i3l(ta)Bl{ta) ' 


^Bl{ta)Bl{h) 


— B\ {ta)B'l[t 


a)- 


- i3l{h)i3\{h) - 


^ Bl{h)Bj{ta))Y 



m W 
2^V 15 



(120) 



Notice that under transformation {xi^X2) — ^ {x^y) both 
W and -D do not change. 



APPENDIX B 

We prove here that both Sd [x] and Xc/ {i) are independent 
of the choice of the fundamental system. The proof can 
be done in two steps. Firstly we show that both S'c/[x] 
and Xc; do not depend on the scaling of U; and . Indeed, 
if we rescale, for example ui — > aui , then 



Ua CtUa 

D3 aDs 



D2 aD2 , 
Z?4 aDi . 



Analogous relations are valid for other vectors. It is now 
simple to see that both Sd and Xc/ remain unchanged 
under such a rescaling. 



Secondly, we show that Sd and Xc/ remain unchanged 
under the substitution 



ui ^ ui + au2 + /?vi + 7V2 , 



(121) 



where a, /?, 7 are arbitrary real constants. The previous 
substitution is possible to achieve, for example, in suc- 
cessive steps: 



Ul 

Vl 

V2 



Ul + /3vi , 



V2 



7 



-U2 



It may be directly seen that both Sd and Xd are invari- 
ant under each of the former substitutions and so they 
are invariant with respect to ( |l2l| ), too. To see it more 
explicitly let us perform, for instance, the substitution 
Ul — > Ul + /3vi, then 



Ua^Ua 

D4 -> D4 



Vb- 
D3 



Plugging the previous substitution into (16), we get that 
^ciit) ~^ Xc/(i). The same is true for Sd as it might be 
directly seen from relations ( p^ and ( |119| ). 

The previous two observations therefore lead us to the 
conclusion that Sd does not depend on the particular 
choice of the fundamental system of solutions (expressible 
as a linear combination of u.^ and v^). 



APPENDIX C 

We derive here some relations used in Section IIIA. 
As we have mentioned in Subsection IIB, having one 
solution, say for (a;i,X2) coordinates, we can get an- 
other one if we multiply the original one by (T3. So 
namely if one has the fundamental system of solutions 
(ui, U2, Vl, V2), one can generate another fundamental 
system (asUi, cr3U2, csvi, cr3V2) (it is simple to see that 
this is indeed a fundamental system by looking at the 
Wronskian). Identical reasonings as in the Subsection 
IID will lead us to the result 



Xc; (t) 



1 

D 



x\ 



Blit) 
Bl{t) 

Blit) 
Blit) 



Xb 



Blit) 
Blit) 

Blit) 
Blit) 



(122) 



Comparing (122) with (|18| ) we get the following useful 
identities: 



Blit) ^ Blit) 
Blit) ^ Blit) 



Blit) = Blit). 
Blit) = Blit). 



(123) 
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Similar analysis may be done with the fundamental 
systems {<T3Ui{-t),a3U2{-t),a3Vi{-t),a3V2{-t)) and 

(l(T2Ui(-t),io-2U2(-t),Z(T2Vi(-t),iO-2V2(-i))- this 

case it can be directly checked that we have 



(124) 



From the definition ( pO[ ) and relations (123) we may ob- 
serve that 



dta 
so namely 



D = 2B\{ta) 



dtb 



D 



dta \ dta 



2 ' ^ —InD = Tr D^^ — D 



D dth 



dth 



(125) 



(126) 



(127) 



Here D — detD. Because D{ta, tb) — D{tb — <q) it fulfills 
the equation 



dD{ta,tb) ^ dD{ta,tb) 
dta dtb 

and so we have the identity 

Blita)^ -Blitb). 

APPENDIX D 



= 0, 



(128) 



We prove here the relation (|3l|). For this purpose it is 
simpler to work in the {x, y) coordinates. The kernel 
can be then constructed in analogous way as in [xi, X2) 
coordinate. A simple calculation shows that 



c r 1 



2 dta 

1— «!) 



i 

2 dtb 



dD 



+ Tan (e"' Blitb) + e"^'"" Blitb)) 
and that 

Blitb)Blitb) = WD. 
The latter allows to identify 



(129) 



Blitb) 
D 



Blitb) 
D 



This identification fixes a modulo iir and leads to the 
equation 



a^iln 



Blitb) 
Blitb), 



(130) 



In addition, from the symmetry reasonings result the fol- 
lowing useful relations 

i3litb)^~i3lita) , i3litb)^-i3lita). 

Let us now consider BUtb). Its explicit structure reads 



Blitb) 



u\itb) u\itb) vlitb) vlitb) 

ujita) ulita) vfita) w|(ia) 
Ui(4) U2(<6) Vi(4) V2(4) 



(131) 



Rules for differentiation of determinants tell that Blitb) 
fulfills the equation 



Blitb) ^0, (132) 



with the boundary condition Blitb)\ta=t 
eral solution of (132) reads 



Blitb 



0. The gen- 
*f^*»/(4)) . (133) 



Here / and / are some functions of tb- Applying the 
boundary condition we get 



Blitb) ^ C e-^'^'^-'^^hinnitb ~ ta) 



(134) 



with C being a constant. The result is clearly the only 
one which is compatible with the differential equation for 
Blita) i^ -Blitb)): 

tiBlita) + ^ ^Blita) + - Blita) = , (135) 

dtj^ m dtb m 

fulfilling the boundary condition -B|(ta)|i[,=ta — 0. 

The same reasonings can be now applied on BUtb). The 
latter fulfills the differential equation 

^ Blitb) + ^ Blitb) + - Blitb) =0, (136) 



dt 



m dta 



with the boundary condition Blitb)\t 
tion is 



0. The solu- 



Blitb) = Ce^f*'-*") sinf7(4 - ta) , 



(137) 



with C being a constant. The result is the only one which 
is compatible with the differential equation for B^ita) 

i^ -Blitb)): 



- 3 4rBlita) + 3 Blita) = , 



m dtb 



(138) 



fulfilling the boundary con di tion B\ita) \tK= t„ = 0. 
Gathering the results (13C), (134) and (137) together we 
obtain 

aitaM) = r(ta - tb) + InC - InC , (139) 
modulo ITT. Using I'Hopital rule one may persuade oneself 
thatlimt^_t, BUtb)/ BUtb) ^ C/C =1. Thus, InC- 
InC is either zero or purely imaginary. 



21 



APPENDIX E 

Let us first emphasize that the formulation of a time re- 
versal transformation must avoid using properties of the 
forces or interactions that determine the dynamics, be- 
cause it is the transformation properties of the dynamic 
equations which we seek to determine. The latter is 
the crux often overlooked by many authors. Since the 
kinematics are those properties of the motion that are 
independent of the dynamics, we require that the "ad- 
missible" time-reversal transformation should be formu- 
lated in kinematic terms. This means that the time- 
reversal transformation must be consistent with the al- 
gebraic structure of the operators representing the (kine- 
matic) observables and that in the absence of forces or 
interactions (i.e., in the absence of causal effects), the 
dynamic equations must be left invariant. 

In our case the kinematic observables may be taken to be 
x" and P" (note that P = mx are the kinetic momenta 
and not the full canonical momenta (pT|)). Working in 
(xi,X2) coordinates, the algebraic structure is then de- 
termined by the Heisenberg-Weyl group 



K,P^]=*n(a3) 



a/3 



[x,a3] - [P,a3] -0, (140) 



On the other hand, if Hq is the Hamiltonian in the ab- 
sence of interaction (7 = 0), the dynamic (Schrodinger) 
equation 



ih-m)) 



Ho\m) , 



must transform under time reversal T into 
d 



(141) 



(142) 



where t' 
obtain 



-t. Applying T to both sides of (141) we 



dt 



(143) 



Comparing (143) with ( |142| ) and using the requirement 
that Hq is invariant under time reversal we obtain the 
relation: TiT~^ = —i. Invoking Wigner's theorem 
the latter implies that T must be antiunitary and thus 
may be written as T = U)C with /C being the complex 
conjugation operator and U being some unitary operator. 

At the same time, in accordance with the classical con- 
ditions, time reversal requires that 



j)X, 



(144) 



The matrix B is supposed to leave Hq invariant under 
the time reversal. This means that B* (0-3)8 — and 
= 1 (i.e., time reversal when repeated must restore 
the original situation). So B must be part of a dis- 
crete (two-element) subgroup of 0(1, 1; R) (i.e., the real 



Lorentz group in the plane). It is well known that the 
only matrices fulfilling the above conditions are ±1 and 
±0-3. 

To decide the form of B we use the fact that Ho is invari- 
ant under 0(1, 1; R). This means that x and P transform 
under 0(1, 1; i?) in the usual manner, i.e.. 



x' = U{e)xU-\e) = G(e)x, 
P' = U{e)PU~\e) =G{e)P , 



(145) 



where U{e) is a unitary representation of 0(1,1;_R) in 
the state space and G(e) is an element of 0(1, l;i?) in 
2-dimensional vector space. As a result, the following 
relation must hold for any G E 0{1,1; R): 



ix. 



(146) 



However, the very same relation may be recast in a 
slightly different form, namely 



= [/BC/"^(Gx : 



(147) 



Comparing both (|l4^) and ( |l47D we get that GB = 
(G*BG^ As this must be true for all G G 0(1,1;P), 
we can choose 



(e) — exp (ecTi) = cosh(e) -I- sinh(e)cri 



(148) 



It is then obvious that the case B = ±1 is ruled out 
and we are left with B = zta^. However, the sign 
is the only plausible one. This is because the signature 
of the time reversal should be preserved under continu- 
ous change of coordinates and so namely when we shrink 
the X2 coordinate into the origin (i.e., perform a dimen- 
sional reduction) xi coordinate must behave under time 
reversal as in ordinary ID l.h.o.. 

As a upshot of the performed analysis we have the fol- 
lowing transformations: 



X - 




= TxT-i 


= 0-3X, 


p 


-> 


= TpT-i 


= - o-sP , 


r — 




= TrT-^ -- 


= r , 


u - 









(149) 



(here p are full, i.e., canonical momenta). Similarly, we 
find that 



TAT-^ = - 


- A, 




= B, 








= -J- 


TJiT-^ = 









(150) 

and the time-reversed commutation relations 
[:r^p^] = -^fi(a3)"^ [x,, 03] = [p„ 03] = . (151) 
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Transformation rules (149) or ( |15C| ) assert that 

rcr-i = c , r JaT-i = Ja . (152) 

We thus finaUy arrive at the conclusion that THT^^ — 
H. The latter is not actually compatible with the time 
reversal presented in Ref. where B = 1 was incor- 
rectly assumed. 

APPENDIX F 

Using the fact that SU{\, 1) ladder operators are J+ and 
J_ and the SU{1, 1) vacuum state is the state | j, (i.e., 
J-^ \ j, I J I) =0), we may write 



I ,s , ,1 V'^'^' / (2|j|)!(m-b1)! 



-2|j| 



J 



(153) 



Here \z)) — exp(zJ-i-) is the (unnormalized) coher- 

ent state of S'C/(1,1) 133] (which can be identified with 
the Gelfand-Neimark z-basis IpOl). 



In deriving Eq.(153) we have employed ( |76| ) together with 
the "annihilation" and "creation" relations 



{J-?\j,m) 



0, 



,fc - k) , m- k> \ j\ 
m — k < \j\ 



{J+)'^\j, TTfi) = bj^yn,k \j, m + k) . 



(154) 



where 





(m 


- j)!(rn 




I 


(m 




— fc)!(rn 


+ J 


-fc)! 


(m 


+ k 


+ j)!(m 




-jV- 




(m 


+ jy.{m 


-j) 





Note that (J+J_) |j,m) = (m^ - j^^ |j,m). Form (|l5^ ) 
suggests that the state |V'^ ;) can be alternatively inter- 
preted as an excited SU{1,1) coherent state |1)). Be- 
cause L" (x ) is a polynomial of n-th order in x, the re- 
lation (153) asserts that there is up to (m — |j|) new 
"SU{1, 1) excitations" condensed into the coherent state 
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